Quantifying the non-Gaussian character of a quantum state by quantum relative entropy 



Marco G. Genoni,''^ Matteo G. A. Paris, ''^-^ and Konrad Banaszek'* 

' Dipartimento di Fisica dell' Universitci di Milano, 1-20133, Milano, Italia 
'CNISM, UdR Milano Universita, 1-20133 Milano, Italia 
^ISI Foundation, 1-10133 Torino, Italia 
Institute of Physics, Nicolaus Copernicus University, PL-87-100 Torun, Poland 
(Dated: January 7, 2009) 

We introduce a novel measure to quantify thie non-Gaussian character of a quantum state: the quantum relative 
entropy between the state under examination and a reference Gaussian state. We analyze in details the properties 
of our measure and illustrate its relationships with relevant quantities in quantum information as the Holevo 
bound and the conditional entropy; in particular a necessary condition for the Gaussian character of a quantum 
channel is also derived. The evolution of non-Gaussianity (nonG) is analyzed for quantum states undergoing 
conditional Gaussification towards twin-beam and de-Gaussification driven by Kerr interaction. Our analysis 
allows to assess nonG as a resource for quantum information and, in turn, to evaluate the performances of 
Gaussification and de-Gaussification protocols. 
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Introduction — The use of Gaussian states and operations 
allows the implementation of relevant quantum information 
protocols including teleportation, dense coding and quantum 
cloning IH]. Indeed, the Gaussian sector of the Hilbert space 
plays a crucial role in quantum information processing with 
continuous variables (C V), especially for what concerns quan- 
tum optical implementations On the other hand, quan- 
tum information protocols required for long distance commu- 
nication, as for example entanglement distillation and entan- 
glement swapping, require nonG operations Besides, it 
has been demonstrated that using nonG states and operations 
teleportation S 01 and cloning |7j] of quantum states may 
be improved. Indeed, de-Gaussification protocols for single- 
mode and two-mode states have been proposed flSSHSl 
and realized ifioll . From a more theoretical point of view, it 
should be noticed that any strongly superadditive and continu- 
ous functional is minimized, at fixed covariance matrix (CM), 
by Gaussian states. This is crucial to prove extremality of 



Gaussian states and Gaussian operations 111 U 11211 for various 



quantities such as channel capacities Ill3ll . multipartite entan- 
glement measures ifTill and distillable secret key in quantum 
key distribution protocols. Overall, nonG appears to be a re- 
source for CV quantum information and a question naturally 
arises on whether a convenient measure to quantify the nonG 
character of a quantum state may be introduced. Notice that 
the notion of nonG already appeared in classical statistics in 
the framework of independent component analysis 115 1. 



The first measure of nonG of a CV state g has been sug- 
gested in lfl6ll based on the Hilbert-Schmidt (HS) distance be- 
tween Q and a reference Gaussian state. In turn, the HS-based 
measure has been used to characterize the role of nonG as a re- 
source for teleportation llVlllSll and in promiscuous quantum 
correlations in CV systems 11911 . Here we introduce a novel 
measure 5[q\ based on the quantum relative entropy between g 
and a reference Gaussian state. The novel quantity is related to 
information measures and allows to assess nonG as a resource 
for quantum information as well as the performances of Gaus- 



sification and de-Gaussification protocols. In the following, 
after introducing its formal definition and showing that it can 
be easily computed for any state, either single-mode or multi- 
mode, we analyze in details the properties of 5[q\ as well as 
its dynamics under Gaussification 12911 and de-Gaussification 
protocols. 

Gaussian states — Let us consider a CV system made of 
d bosonic modes described by the mode operators a^, k = 
\ . . .d, with commutation relations [a/;, aj] = 5k j- A quan- 
tum state g of d bosonic modes is fully described by its char- 
acteristic function x[gi](A) = Tr[g D{\)] where D{\) = 
^k^i Dk{Xk) is the d-mode displacement operator, with A = 
(Ai, . . . , Ad)^, Afc G C, and where Dk{\k) = exp{Afca[, - 
A^afc} is the single-mode displacement operator. The canoni- 
cal operators are given by = {a k + al,) /y/2 and pk — {ak^ 
a\)/\f2i with commutation relations given by {qj,Pk\ = i^jk- 
Upon introducing the vector R = (qi,pi, . . . ,qd,Pd)^, the 
CM cr = (7[g] and the vector of mean values X = X[g] of 
a quantum state g are defined as akj = \{RkRj + RjRk) — 
{Rj){Rk) and = (Rj), where (O) = Tr[g O] is the ex- 
pectation value of the operator O. A quantum state go is said 
to be Gaussian if its characteristic function is Gaussian, that is 
x[^?G](A) = cxp | — iA^^crA + X^a|, where A is the real 

vector A = (ReAi,ImAi, . . . , ReA^, ImAd)^. Once the CM 
and the vectors of mean values are given, a Gaussian state is 
fully determined. For a system of d bosonic modes the most 
general Gaussian state is described by d(2d + 3) independent 
parameters. 

Non-Gaussianity of a CV state — The von Neumann en- 
tropy of a quantum state is defined as S{g) = —TT[g \ogg]. 
The von Neumann entropy is non-negative and equals zero 
iff g is a pure state. In order to quantify the nonG char- 
acter of a quantum state g we employ the quantum relative 
entropy (QRE) S(f)||T) = Ti[g{\ogg — logr)] between g 
and a reference Gaussian state r. As for its classical coun- 
terpart, the Kullback-Leibler divergence, it can be demon- 
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strated that < S(p||t) < oo when it is definite, i.e. when 
supp g C supp r. In particular S(f3||T) = iff g = t. This 
quantity, though not defining a proper metric in the Hilbert 
space, has been widely used in different fields of quantum 
information as a measure of statistical distinguishability for 
quantum states |23, 211- Therefore, given a quantum state g 
with finite first and second moments, we define its nonG as 
S[g] = S(£i||t), where the reference state r is is the Gaussian 
state with X[g] = X[t] and cr[g] = (t[t], i.e. the Gaussian 
state with the same CM a and the same vector X of the state 
g. Finally, since r is Gaussian, then log r is a polynomial op- 
erator of the second order in the canonical variables which, 
together with the fact that r and p have the same CM leads to 
Tr[(T - g) logr] = 10, i.e. S(g||r) = S(r) - S{g). Thus 
we have 



S[g] = S(r) - Sig) 



(1) 



i.e. nonG is the difference between the von Neumann en- 
tropies of r and g. In turn, several properties of the nonG 
measure Mo] may be derived from the fundamental properties 
of QRE i20l l2in . In the following we summarize the relevant 



ones by the following Lemmas: 

LI: 6[g] is a well defined non negative quantity, that is 
< ^[p] < oo and S[g] ~ Q iff p is a Gaussian state. Proof: 
Nonnegativity is guaranteed by the nonnegativity of the quan- 
tum relative entropy. Moreover, if (5[f3] = then g ~ t and 
thus it is a Gaussian state. If p is a Gaussian state, then it is 
uniquely identified by its first and second moments and thus 
the reference Gaussian state r is given hy t = g, which, in 
turn, leads to 5[g\ — S{g\\T) = 0. 

L2: 5[g\ is a continuous functional. Proof: It follows from the 
continuity of trace operation and QRE. 
L3: 5[g\ is additive for factorized states: 5[gi (g) ^2] = 5[gi\ + 
5[g2\- As a corollary we have that if g2 is a Gaussian state, 
then 5 [^?] = S[gi\. Proof: The overall reference Gaussian state 
is the tensor product of the relative reference Gaussian states, 
T = Ti ® T2. The lemma thus follows from the additivity of 
QRE and the corollary from LI. 

L4: For a set of states {gu} having the same first and sec- 
ond moments, then nonG is a convex functional, that is 

'^EfcPfcfi'fc] < Y.kPk6[gi], with Y^kPk = 1- Proof: The 
states gk, having the same first and second moments, have 
the same reference Gaussian state t which in turn is the 
reference Gaussian state of the convex combination g = 
TlikPkgk- Since conditional entropy S(£i||r) is a jointly con- 
vex functional respect to both states, we have ^^[^i^Pkgk] = 
S(EfcPfc£'fe||'r) < Y.kPkS{gk\\T) = J2kPk^l^k]- □ Notice 
that, in general, nonG is not convex, as it may easily proved 
upon considering the convex combination of two Gaussian 
states with different parameters. 

L5: If Ub is a unitary evolution corresponding to a symplectic 
transformation in the phase space, i.e. if Ub = cxp{—iH} 
with H at most bilinear in the field operator, then S[UbgUl] = 
5[g]. Proof: Let us consider g' = UbgU^, where U is at 
most bilinear in the field-mode, then its CM transforms as 



(T[g'] = T,(T[g]'S'^ , S being the symplectic transformation as- 
sociated to U. At the same time the vector of mean values 
simply translates to X' = X + Xq. Since any Gaussian state 
is fully characterized by its first and second moments, then the 
reference state must necessarily transform as r' = UbTl/^, i.e. 
with the same unitary transformation U. Lemma follows from 
invariance of QRE under unitary transformations. 
L6: nonG is monotonically decreasing under partial trace, that 
is i5[TrB[£i]] < 6[g]. Proof: Let us consider g' — Tislg]- Its 
CM is the submatrix of cr[g] and its first moment vector is 
the subvector of X[g\ corresponding to the relevant Hilbert 
space. As before, also the new reference Gaussian state must 
necessarily transform as r' = TrB[r]. QRE is monotonous 
decreasing under partial trace and the lemma is proved. 

L7: nonG is monotonically decreasing under Gaussian quan- 
tum channels, that is < 5[g\. Proof: Any Gaus- 
sian quantum channel can be written as £g(p) = Tr^ [J7f,(g (g) 
te)uI], where Ub is a unitary operation corresponding to an 
Hamiltonian at most bilinear in the field modes and where te 
is a Gaussian state ll23ll . Then, by using lemmas L3, L5 and 
L6 we 

obtain 5 [£G(e)] < 5[Ub{g®TE)Ul] ^ 5[g\. □ In turn, 
L7 provides a necessary condition for a channel to be Gaus- 
sian: given a quantum channel £, and a generic quantum state 
g, if the inequality (5[£(gi)] < 5[g\ is not fulfilled, the channel 
is nonG. 

Maximally non-Gaussian states — Let us now consider a 
single mode [d = 1) system and look for states with the max- 
imum amount of nonG at fixed average number of photons 



N 



Since 5[g] = S(t) — S(g) we have to maxi- 
mize S(r) and, at the same time, minimize S{g). For a single- 
mode system the most general Gaussian state can be written as 
ga ^ D{a)S{C)v{nt)S^C)DHa), D{a) being the displace- 
ment operator, S{Q = exTp{^(a^'^ — ^(*a^) the squeezing 

operator, e C, and i/(nt) = (1 + nt)^^[nt/{l + nt)]"^" 
a thermal state with nt average number of photons. Dis- 
placement and squeezing applied to thermal states increase 
the overall energy, while entropy is an increasing monotonous 
function of the number of thermal photons nt and is invariant 
under unitary operations, thus, at fixed energy, S(r) is max- 
imized for r — v{N). Therefore, the state with the maxi- 
mum amount of nonG must be a pure state (in order to have 
S{g) = 0) with the same CM cr = (A^ + i)I of the ther- 
mal state v{N). These properties individuate the superposi- 
tions of Fock states \iIjn) ~ '^k^k\n + Ik) where n > 0, 
^k > Ik-i + 3 or = 0, with the constraint N = {a^a), 
i.e n + Y.k W?h = N = {dct (T[iy{N)])i - i. These 
represent maximally nonG states, and include Fock states 
IV'jv) = |iV) as a special case. Let us consider now d- 
mode quantum states with fixed average number of photons 
X]fe=i Tr[aj^afcp] = N = iik- Also in this case maximally 
nonG states are pure states; the CM being equal to that of a 
multimode classical state t = R ®k v{mk)R} , "^/^ ~ 
where we denote with R a generic set of symplectic passive 
operations {e.g. beam splitter evolution) which do not increase 
the energy. In order to maximize S(r) = Tlik-'i^i'^k)) we 
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have to choose mfc = N/d for every k. As for example, fac- 
torized states of the form \9 n) ~ iV'JV/d)®''- whose reference 
Gaussian states are r ~ [v{N/d)]®^'-, are maximally nonG 
states at fixed N . Of course for the multi-mode case there are 
other more complicated classes of maximally nonG states that 
include also entangled pure states. Finally, we observe that the 
maximum value of nonG is a monotonous increasing function 
of the number of photons N . 

Non-Gaussianity in quantum information — Gaussian states 
are extremal for several functionals in quantum information 
1 1 in . In the following we consider two relevant examples, 
and show how extremality properties may be quantified by the 
nonG measure 6[g\. Let us first consider a generic communi- 
cation channel where the letters from an alphabet are encoded 
onto a set of quantum states Qk with probabilities pk- The 
Holevo Bound represents the upper bound to the accessible 
information, and is defined as xio) — S(f?) — J2kPk S(^?fc) 
where g = J^kPkQk is the overall ensemble sent through 
the channel. Upon fixing the CM (and the first moments) 
of g we rewrite the Holevo bound as x{q) ~ S(t) — 6[g] — 
J2kPk S(gfe), where t is the Gaussian reference of g. This 
highlights the role of the nonG d[g] of the overall state in de- 
termining the amount of accessible information: at fixed CM 
the most convenient encoding corresponds to a set of pure 
states gk, S{gk) ~ 0, forming an overall Gaussian ensem- 
ble with the largest entropy. In other words, at fixed CM, 
we achieve the maximum value of x upon encoding encoding 
symbols onto the eigenstates of the corresponding Gaussian 
state 12411 . If the alphabet is encoded onto the eigenstates of a 
given state g, we have x{q) = S(t) — S[g]. This suggests an 
operational interpretation of nonG S[g] as the loss of informa- 
tion we get by encoding symbols on the eigenstates of g rather 
than on those of its reference Gaussian state. 

Let us now consider the state gAB describing two quan- 
tum systems A and B and define the conditional entropy 
S{A\B) = S{gAB) - Sign)- Let us fix the CM of gAB 
and thus also that of gB, and consider the reference Gaus- 
sian states TAB and tb- We may write S{A\B) ~ Sg{A\B) — 
{S[gAB] - S[gB]) where Sg{A\B) = S{tab) - S{tb), i.e. 
the conditional entropy evaluated for the reference Gaussian 
states TAB and tb- Then, upon using L6 we have (5[gA_B] — 
6[gB] > and thus S{A\B) < Sg{A\B), i.e the maximum 
of conditional entropy at fixed CM is achieved by Gaussian 
states. In classical information theory the conditional entropy 
H{X\Y) = if (X, y)— 77(F), where von Neumann entropies 
are replaced by Shannon entropies of classical probability dis- 
tributions, is a positive quantity and may be interpreted lEsll 
as the amount of partial information that Alice must send to 
Bob so that he gains full knowledge of X given his previous 
knowledge from Y. When quantum systems are involved the 
conditional entropy may be negative, negativity being a suffi- 
cient condition for the entanglement of the overall state gAB ■ 
This negative information may be seen as follows lE6ll for a 
discrete variable quantum system. Given an unknown quan- 
tum state distributed over two systems, we can discriminate 
between two different cases: ifS(j4|_B) > 0, as in the classical 



case, it gives the amount of information that Alice should send 
to Bob to give him the full knowledge of the overall state gAB. 
When S{A\B) < Alice does not need to send any informa- 
tion to Bob and moreover they gain —S{A\B) ebits. If we 
conjecture that this interpretation can be extended to the CV 
case, the relation S(A|i?) < Sg{A\B) ensures that, at fixed 
CM, nonG states always perform better: Alice needs to send 
less information, or, for negative values of the conditional en- 
tropy, more entanglement is gained. Moreover, since nega- 
tivity of conditional entropy is a sufficient condition for en- 
tanglement II27II we have that for any given bipartite quantum 
state gAB, if the conditional entropy of the reference Gaussian 
state TAB is negative, then gAB is an entangled state. Though 
being a weaker condition than the negativity of S{A\B), this 
is a simple and easy computable test for entanglement which 
is equivalent to evaluate the symplectic eigenvalues ll28ll of the 
involved Gaussian states. 

Gaussification and de-Gaussification protocols — Since the 
amount of nonG of a state affects its performances in quantum 
information protocols a question naturally arises on whether 
this may be engineered or modified at will. As concerns Gaus- 
sification, Lemma L7 assures that Gaussian maps do not in- 
crease nonG. In turn, the simplest example of Gaussification 
map is provided by dissipation in a thermal bath iflill . which 
follows from bilinear interactions between the systems under 
investigation and the environment. On the other hand, a condi- 
tional iterative Gaussification protocol has been recently pro- 
posed ligll which cannnot be reduced to a trace-preserving 
Gaussian quantum map. It requires only the use of passive el- 
ements and on/off photodetectors. Given a bipartite pure state 
in the Schmidt form j^/jC'')) — Yl^=a'^rwn\ri,n) the state at 
k+ 1-th step of the protocol has the same Schmidt form with 



)arlrO:"j^'_r n-r- have Considered 



(fe+l) _ „_„ 
Ltn.n — ^ / n 

the initial nonG superposition = (1 + A^)"^/2(|0, 0) + 

A|l, 1)) which is asymptotically driven towards the Gaussian 
twin-beam state |?/') = Vl — A^ J2^=o -^"I^^j "-)■ We have 
evaluated nonG at any step of the protocol, for every value of 
A. 
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FIG. 1: (Left): nonG after some steps of the conditional Gaussifi- 
cation protocol of Ref. i29il considering as the initial state the nonG 
superposition l?/)'"^) = (1 + A^)"^''^(10, 0) + A|l, 1>): black-solid: 
initial state; black-dashed : step 1; gray-solid: step 2; gray-dashed: 
step 3. (Right): black-solid: initial state; black-dashed : step 5; gray- 
solid : step 10; gray-dashed: step 20. 

Results are reported in Fig. [T] for the first steps, nonG de- 
creases monotonically for almost all values of A (only at the 
third step, for Awl the state is more nonG than at the pre- 
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vious steps). Notice that increasing the number of steps nonG 
may also increase, e.g, for A w 1, (5 reaches very high val- 
ues and the maximum value increases. On the other hand, the 
overall effectiveness of the protocol is confirmed by our anal- 
ysis, since the range of values of A for which (5 « increases 
at each step of the protocol. In other words, though not be- 
ing a proper Gaussian map, the conditional protocol of ll29ll 
indeed provides an effective Gaussification procedure. 

Conditional de-Gaussification procedures have been re- 
cently proposed and demonstrated lUSIllESl- Here we rather 
consider the unitary de-Gaussification evoluti on p rovided by 
self-Kerr interaction Uj = cxp{— 17(0^0)^} BSCj Bill , which 
does not correspond to a symplectic transformation and leads 
to a nonG state even if applied to a Gaussian state. We have 
evaluated the nonG of the state obtained from a coherent state 
I a) undergoing Kerr interaction. Results are reported in Fig. |2] 
as a function of the average number of photons (up to 10^ pho- 
tons) and for different values of the coupling constant 7 . As it 
is apparent from the plot, nonG is an increasing function of the 
number of photons and the Kerr coupling 7. For 7 « 10^^, 
the maximum nonG achievable at fixed energy is quite rapidly 
achieved. For more realistic values of the nonlinear coupling, 
i.e. 7 < 10^^ nonG states may be obtained only for a large 
average number of photons in the output state. In fact, to ob- 
tain entanglement, experimental realizations 1 3^ 31 ] involve 
pulses with an average number of the order of 10^ photons, 
which are needed to compensate the almost vanishing small 
Kerr nonlinearities of standard glass fiber. We finally no- 
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FIG. 2: NonG of coherent states undergoing Kerr interaction as a 
function of the average number of photons and for different values 
of the coupling constant 7. Black dashed line: 7 = 10^^; black 
dotted: 7 = 10"'*; solid black: 7 = 10^'^. The gray solid lines is 
the maximum nonG at fixed number of photons. 

tice that a good measure for the nonG character of quantum 
states allows us to define a measure of the nonG character of 
a quantum operations. Let us denote by Q the whole set of 
Gaussian states. A convenient definition for the nonG of of a 
map £ reads as 6[£] = maxg^g S[£{g)], where £{g) denotes 
the quantum state obtained after the evolution imposed by the 
map. 

Conclusions — We have introduced a novel measure to 



quantify the nonG character of a CV quantum state based on 
quantum relative entropy. We have analyzed in details the 
properties owned by this measure and its relation with some 
relevant quantities in quantum information. In particular, a 
necessary condition for the Gaussian character of a quantum 
channel and a sufficient condition for entanglement of bipar- 
tite quantum states can be derived. Our measure is easily 
computable for any CV state and allows to assess nonG as 
a resource for quantum technology. In turn, we exploited our 
measure to evaluate the performances of conditional Gaussi- 
fication towards twin-beam and de-Gaussification processes 
driven by Kerr interaction. 
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